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Abstract 

Hamiltonian systems are one of the most important class of dynamical systems with a geometric 
structure called symplecticity and the numerical algorithms which can preserve such geometric 
structure are of interest. In this article we study the construction of symplectic (partitioned) 
Runge-Kutta methods with continuous stage, which provides a new and simple way to construct 
symplectic (partitioned) Runge-Kutta methods in classical sense. This line of construction of 
symplectic methods relies heavily on the expansion of orthogonal polynomials and the simplifying 
assumptions for (partitioned) Runge-Kutta type methods. 
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1. Introduction 


Geometric numerical integration is a subfield of the numerical solution of differential equations, 
and it turn out to be very efficient for long-time simulating the dynamical behavior of those systems 
with special structures |14l |. For Hamiltonian systems, it is known that symplecticity of the flow is 
an important characteristic property of such systems [l| and to correctly reproduce this property, 
the numerical algorithms called symplectic methods have been proposed and extensively studied 
(see B 0, Q and references therein), among of which symplectic (par titioned) Runge-Kutta 
methods are one of the most important subclass of such algorithms 20|, 124,122, |23( . In this article, 
we focus on a new and simple way to construct these methods. 

Runge-Kutta (RK) method with continuous stage has been discussed and investigated by sev¬ 


rn 


eral authors recently |l5l. l25l. l26l. I28l.ll8l|. the idea of which was firstly presented by Butcher in 1970s 
There are few literatures involving this kind of methods since Butcher’s seminal work, until 
recently, [l5| exploits it to interpret the energy-preserving collocation methods. Subsequently, [25 1 
investigates four kinds of time finite element methods and relates them to Runge-Kutta methods 


with continuous stage. Moreover, it is shown in 251] that some energy-preserving_RK methods 
including s-stage trapezoidal method 16], Hamiltonian boundary value methods B B, average 
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vector field method [l9|, can also be related to such methods. [26|, |28| discuss the construc¬ 
tion of Runge-Kutta methods with continuous stage and present the characterizations for sev¬ 
eral geometric-structure preserving methods including symplectic methods, symmetric methods 
and energy-preserving methods. Exponentially-fitted continuous-stage R ung e-Kutta methods with 
energy-preserving property for Hamiltonian systems are also presented in 1^. In this paper, we aim 
at further investigation of construction of symplectic (partitioned) Runge-Kutta methods, along 
the line of expansion of orthogonal polynomials and with the help of simplifying assumptions for 
(partitioned) Runge-Kutta type methods. This approach can naturally provides a simple way to 
construct symplectic (partitioned) Runge-Kutta methods in classical sense. 

The outline of the paper is as follows. In Section 2, we provide some preliminaries for (par¬ 
titioned) Runge-Kutta methods with continuous stage and review some existing results for the 
construction of general RK type methods with continuous stage. In section 3, we discuss the con¬ 
struction of symplectic (partitioned) Runge-Kutta methods with continuous stage as well as the 
symplectic (partitioned) Runge-Kutta methods in classical sense. At last, we conclude this paper. 


2. Construction of (P)RK type methods 

In this section, we first introduce the definition of (partitioned) Runge-Kutta method with 
continuous stage by following the formulation proposed in 1^, then show some results which are 
useful in constructing (P)RK type methods. 

2.1. Continuous-stage RK method 

Consider a first-order system of ordinary differential equations (ODEs) 


( 2 . 1 ) 

is a sufficiently 


fz = /(t,z), 

\z{to) = zo G 

where the upper dot indicates differentiation with respect to t and / : M x —? 

smooth vector function. 

Definition 2.1. [l^. Let Ar^a be a function of two variables t, a ^ [0,1], be a function of 
r G [0,1]. Define At^^t da. The one-step method <l>h : zq zi given by 

Zr = Zo-hh [ Ar,af{to + C^h, Z„) d(T, T G [0, 1], 

Jo 

Zi = Zo -b h / Brf{to + Crh, Zr) dr, 

Jo 

is called a continuous-stage Runge-Kutta (csRK) method, where Zr ~ z{to + Crh), zi Ri z(to + h). 


( 2 . 2 ) 


In this paper, the construction of csRK methods mainly relies on the following simplifying 
assumptions of order conditions 


1 


^(0: / BrCr'dT = -, K = l,...,^, 

Jo ^ 


1 


Ar,.C^-^da = -C^, K = l,...,rj, 
lo ^ 


€{7]) : / 

Jo 

C 1 

D(C): / BrC^-^Ar,„dT = -B„{l-C:), k = 

Jo ^ 
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Analogously to the classical result given by Butcher in 1964, we can determine the order of a 
csRK method with the help of these simplifying assumptions. 

Theorem 2.1. If a csRK method (j2.2l) with coefficients {At, a, Bt, Ct) satisfies ^{p), ^{a) and 
D{I3), then it is at least of order 


min(/9, 2a + 2, a + /3 + 1). 

Proof. This statement can be proved by using the similar idea to the one used in Theorem 7.4 

3 . □ 


According to Proposition 2.1 shown in 2^, we assume Bt = 1 in the current paper and to 
proceed conveniently, hereafter we also assume Ct = t. In such a case, 1B(|) is reduced to 


f\.-. 

Jo 


dr=-, K = 1,...,C 


K 


which implies holds with ^ = 00 . Now we introduce the i-degree normalized shifted Legendre 
polynomial denoted by Pi{x), which can be explicitly computed by the Rodrigues formula 


Po{x) = 1, P,{x) = 


^/2^TT C 

l\ dx'- 


(x^x-ir), 6 = 1 ,2,3, 


A well-known property of Legendre polynomials is that they are orthogonal to each other with 
respect to the inner product on [0, 1] 


/ P,{t)PK{t)dt = 6 , k = 0 , 1 , 2 , • • • , 

40 

and they as well satisfy the following integration formulae 

[ Po{t) dt = ^iPi{x) + ]:Po{x), 

Jo ^ 

[ P,(t)dt = ^,+iP,+i(x)-^,p,_i(x), 6 = 1,2,3, 

4o 


(2.3) 


[ P,{t)dt = 6,o- f P,{t)dt, 6 = 0,1,2, 
Jx Jo 


where ^ and 6ij is the Kronecker symbol. 

With the help of (|2.3jl . by exploiting the expansion of Legendre polynomials for the coefficients 
of a csRK method and the corresponding simplif ying assumptions, we have the following result for 
construction of csRK methods (more details see [23, [2^ ). 


Lemma 2.1. 


'2711 For a csRK method with R,- = 1 and Ct = t, we have the following statements: 


(a) €{a) is equivalent to the fact that At, a takes the form 


Oi—l 


= PL{x)dxP,{a)+ '^J,{T)P,{a) 


i=0 


where 7 t(r), 6 > a are arbitrary functions; 
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(b) D(/3) is equivalent to the fact that o- takes the form 


A.r„ = 



where At((T), i> (5 are arbitrary functions. 


Theorem 2.2. 128 1 For a csRK method with = 1 and Cr = t (then *B(oo) holds), the following 
two statements are equivalent: 

(a) Both ^{a) andT){f3) hold; 

(b) The coefficient a has the form 

Ni 

-I- 

2 


Ni N2 

(2.4) 


L=0 


L=0 


i>l3 

j>OL 


where Ni = max(a — 1, /3 — 2), N 2 = max(a — 2, /3 — 1), ^ and jij are arbitrary 

real numbers. 

By combining Theorem 12.II and Theorem 12.21 we can easily construct lots of csRK methods, the 
order of which are given by min(oo, 2a + 2, a + /3 + 1) = min(2a + 2, a + j3 + 1). 


2.2. Continuous-stage PRK method 
For a partitioned ODEs given by 


y = f{t, y, z), y{to) = z/o G 

i = g{t, y, z), z{to) = zq G W, 


(2.5) 


we can consider methods that deal with the y and z variables by two different csRK methods, 
which then produce the following definition. 


Definition 2.2. 

one-step method 


271 . 1281] Assume that Cr 



da and Cr 


Jq Ar^a da. The following 


Yr = yo + fi 

Zr = Zq -\- h 


yi = yo + h 
zi = ZQ + h 


[ Ar,crfito + Ccrh,Ycr,Z^)da, rG[0,1], 

Jo 

/ Ar,a9{to + Cah, Y„, Z„) dcT, r G [0,1], 

Jo 

[ Brf{to + Crh,Yr,Zr)dT, 

Jo 


Brg{to -\- Crh, Yr, Zr) dr. 


is referred to as continuous-stage partitioned Runge-Kutta (csPRK) method. 


( 2 . 6 ) 
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Similarly as those order conditions considered for classical PRK method [8(], we propose the 
following simplifying assumptions of order conditions for the csPRK method 


BiC): / BrC^-^CrdT= -, l<K + t<e; 

Jo n + t 


Civ)-- f 


^K-\-L 


Ar^crC^ da =1<K + L<r], tG[0, 1]; 


C{v): 


■ n+L 


'0 


K + i 


I < K + i <f], r G [0,1]; 


R ll _ 

V{C): / BrC^-^C'r A,^^dT = ^ -^ l<K + i<Q, aG[0,l]; 


K + i 


^(C): f 

Jo 


BrC^-^Cr Ar,adT = 


BM - a 


- + 


K + L 


1 < K + L < C, (T G [0, 1]. 


(2.7) 


Theorem 2.3. If a csPRK method (j2.6h with coefficients (At, o-, Br, Ct] Ar,a, Br, Cr) satisfies 
Br = Br for T G [0,1], and moreover, all B{p), C(a), C{a), 'D{/3), V{fd) hold, then it is at least of 
order 

min(/9, 2a + 2, a + /3 + 1). 


Proof. The proof is the same as that for classical (non-partitioned) Runge-Kutta methods (see 
12l|). in which the bi-colored trees have to be considered instead of trees. □ 


If we let Br = Br ^ I and Cr = Cr = r, then by using Theorem 12.21 it is easy to design At, a 
and At,(t such that all the conditions described in Theorem 12.31 are fulfilled. This allows us to 
construct lots of csPRK methods of arbitrarily high order. 


2.3. Classical (P)RK method retrieved 

Though we have discussed the construction of csRK and csPRK methods above, it should 
be noticed that in general the practical implementation of cs(P)RK methods needs the use of 
numerical quadrature formula. Therefore, next we present classical (P)RK methods retrieved by 
using quadrature formulae. 


2.3.1. RK method retrieved 

Applying a quadrature formula denoted by (6j, to (12.2p . we derive an r-stage classical RK 

method 


Zi = zo + h'^{bjAci,cj)fito + Cjh,Zj), i = !,■■■ ,r, 
i=i 

r 

Zi = Zo + h'^ bif{to + Cih, Zi), 
i=l 


( 2 . 8 ) 


where Ac^^cj = {Ar^a)\T=ci,a=cj and Zi « Zc- = (ZT)|T=ci- The scheme (|2.8I) can also be formulated 
by the following Butcher tableau 
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c 

Cj )'rxr 


bT 


where c = (ci, • • • , Cr)'^, 6 = (6i, • • • , br)'^. 


Theorem 2.4. Assume Ar^a is a bivariate polynomial of degree dA in r and degree d^ in a, and the 
quadrature formula is of orde^p. If a csRK method (|2.2I1 with coefficients {Ar^a, Bt, Cr) 

satisfies = 1, Cr = r (then 5S(oo) holds) and both ^{rj), S(C) hold, then the classical RK method 
with coefficients {bjAc^^cj,bi, Ci) is at least of order 


mm(p, 2a + 2,a + j3 + 1), 


where a = mm(? 7 ,p — d”") and (3 = min(C,p — dA). 

Proof. Since fg f{x) dx = ^ bif{ci) holds exactly for any polynomial f{x) of degree up to p — 1, 

i=l 

using the quadrature formula {bi,Ci)'(^i to compute the integrals of Ci(r/), 2)(C) it gives 


1 

'^bicf~^ = -, K = ,p, 

K 


i=l 


'^{bjAci,cj)Cj ^ i = !,■■■ ,r, k = !,■■■ ,a, 

i=i 

^bicf~\bjAc^,c,i) = ^{l-c')), j = !,■■■ ,r, ,(3. 

i=l 


where a = min(7/,p — d^) and j3 = min(^,p — d'^). Note that in the last formula, we have multiplied 
a factor bj afterwards from both sides of the original identity. These formulae imply that the RK 
method with coefficients ^ satisfies B{p), C{a) and D{/3) (i.e., the classical simplifying 

assumptions for RK method, see [a^llSl), which provides the order of the method from a classical 
result by Butcher B- □ 


^The numerical quadrature formula denoted by (bi, Ci)r=i is called order p if and only if f) f{x) da: = ^ bif{ci) 

i=l 

exactly holds for any polynomial f{x) of degree up to p — 1. 
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2.3.2. PRK method retrieved 

Assume that Br = Bt- = I and Cr = Cr = t, applying a quadrature formula denoted by 
(6j,Ci)[^i to (j2.6p . we then derive a r-stage classical PRK method 


Yi = yo + h^{bjAc,,cj)f{to + Cjh,Yj,Zj), i = !,■■■ ,r, 

r 

Zi = zo + h'^(bjXi,cj)g(to + Cjh,Yj,Zj), i = !,■■■ ,r, 
j = ^ 

r 

yi = yo + h'^bif{to + Cih,Yi, Zi), 

i=l 

r 

zi = ZQ + h'^big{to +Cih,Yi,Zi), 
i=l 


(2.9) 


which can be formulated by a pair of Butcher tableaux below 


c 

{bjAc^^Cj)rxr C 

(Pj-^Ci, Cj )r xr 


b^ 



where c = (ci, • • • , Cr)'^, b = {bi, - ■ ■ , br)"^. 

Theorem 2.5. Assume A^^a is a bivariate polynomial of degree dT^ in r and degree in a, 
Ar^a is a bivariate polynomial of degree d'^ in r and degree d^ in a, and the quadrature formula 
ih,Ci)i^i is of order p. If Br, Cr', Ar^a, Br, Cr) satisfies Br = Br ^ 1, Cr = Cr = t (then 

B(oo) holds) and allC{a), C{a), V{fi), B{fi) hold, then the classical PRK method with coefficients 
{bjAc^^cj,bi, Ci', bjAc^^cj,bi, Ci) is of order 

min(p, 2ri + 2,r] + ( + 1), 

where p = min(Q!,p — d^,p — d^) and t) = min(/3,p — d'^,p — d'^). 

Proof. The proof of the current theorem is similar to that of Theorem 12.41 and one needs using a 
known result presented in Q| (see Theorem 4.1, Page 97). □ 


3. Construction of symplectic (P)RK methods 

We consider Hamiltonian differential equations written in the following form 


p = -VqU{p,q), 

q = VpH(p,q), 


(3.1) 


where the Hamiltonian namely the total energy Yi{p,q) is a smooth scalar function of position 
coordinates p G and momenta q G and we specify the initial condition as {p{to), q{to)) = 
{po, Qo) £ It is well-known that the energy of the system (13. ip is a conserved quantity along 
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the analytic solution of the system, and the system has a geometric structure called symplecticity 
Q, IS) 10; 11], i-e., the flow tpt of the system satisfies 


= J, 


where J is a canonical structure matrix 1^ and denotes the derivative of pt with respect to 
initial values. An important numerical integrator for this special system is the so-called symplectic 
method which has an excellent long-time behavior in simulating the exact flow of the 

system. A one-step method '■ (Po, Qo) (Pi, Qi) for solving Hamiltonian systems is called 
symplectic if 

JcD; = J. 

More surveys and discussions about this special class of algorithms can be found in |n|,Q,|2il and 
references therein. 

It is known that for an r-stage RK method with coefficients when applied to the 

Hamiltonian system m, it can preserve the symplecticity of the system only if the RK coefficients 
satisfy 20, 17, 241 

haij + bjUji = bibj, i,j = !,■■■ ,r, (3.2) 

and for an r- stag e PRK method with coefficients {aij,bi,Ci;aij,bi,Ci), the conditions for symplec¬ 
ticity are 23, 14] 


biCLij bjOiji — b^bj, ijj — 1, ■ ■ ■ , r, 
bi = bi, i = !,■■■ ,r. 


(3.3) 


Analogously to the classical results above, for a csRK method with Br = 1, Cr = t, it is symplectic 
if 

(T + Ao-, r = 1, £ [0, 1], (3.4) 


while for a csPRK method with R,- = R,- = 1, C,- = Cr = r, it is symplectic if 

Ar, (T + Act, r = 1, 'T, <7 € [0, 1]. 


(3.5) 


Next, let us proceed by considering the constrnction of symplectic RK and symplectic PRK 
type schemes, separately for both non-classical and classical cases. 


3.1. Symplectic RK type method 

Theorem 3.1. The csRK method (12.2j) with Br = 1, Ct = t is symplectic if Ar,a has the form 

Act,ct = ^+ ^ aijPi{T)Pj{a), Oij G M, (3.6) 

o<i+jei, 

where aij is skew-symmetric, i.e., aij = —otji, i -|-j > 0. 

Proof. Assume A,-,o- can be expanded as a series in terms of basis {PiiT)Pj{cr)}P^Q, written in the 
form 

Ar,a= ^ aijPi{T)Pj{a), Oij G M, 
o<i,jez 
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and thus we have 


Aa,T= ^ aijPi{a)Pj{T) = ^ ajiPj{a)Pi{T). 

Substituting the above two expressions into ()3.4p and collecting the like basis, it gives 

0^00 — 2’ * ~b j ^0, 

which completes the proof. □ 

It is evident that combing Theorem 13.11 with Theorem 12.21 we can construct symplectic csRK 
methods of arbitrary order. Moreover, when consider using a quadrature formula, we have the 
following corollary. 

Corollary 3.1. If the coefficient Ar, o- of a csRK method with Br = 1, Cr = t are given by (13.611 . 
then the RK scheme {bjAc^^cj,bi, ci) obtained by using a quadrature formula {hi,Ci)\^i is always 
symplectic and the order of the method can be determined by Theorem fS.fl 

Proof. When applying a quadrature formula ibi,Ci)\^i to calculate the integrals of the associated 
csRK method, it can be seen that 

^Ci, cj + = 1) b J = 1) ■ ■ ■ ) 

thus 

bi{bjAc^^cj) + bjibiAcj^a) — bibj, i,j = 1, • • • ,r, 

which implies that the conditions for symplecticity, namely (13.211 . are fulfilled. □ 

Example 3.1. Let = 1, C,- = t and take 

-P,+i(u)P,(r)), (3.7) 

L=0 

by using some suitable quadrature formulae, we can retrieve a lot of classical symplectic RK methods 
including {s + l)-stage Lobatto HIE, Radau IB, Radau IIB and Gauss-Legendre RK schemes, and 
the orders are 2s, 2s + 1, 2s + 1, 2s + 2 respectively. 

If we introduce a parameter A by adding a term X^s+i{Ps+i{T)Ps{or) — Ps+i(c’')Ps('r)) to (|3.7n . 
and consider s = 1, then by using 3-point Gaussian quadrature formula it produces a new class of 
symplectic RK methods, which is shown in Table WH 

3.2. Symplectic PRK type method 

For convenience, we first give the following lemma. 

Lemma 3.1. Assume B^ = 1 and = t. If (^r ,cd Bt-, Ct) satisfies £(ry) and T>(C), define 
Ar^a := 1 — Aa^r, then {At, a, Br, Cr) Satisfies £(C) and D{rj). 

Proof. The proof can be directly deduced from Theorem 12.21 □ 
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VIE 


5 

2 

(2+A)^/T5 

5 

(5-2A)\/T5 

10 


36 

9 

45 

36 

90 


^ + 
36 ^ 

(2+A)^/T5 


2 

5 

(2+A)^/T5 

) 

72 


9 

36 

72 

VTE 

^ + 
36 ^ 

{5-2\)VP5 

^ + 
9 ^ 

(2+A)^/T5 


5 

10 

90 

45 


36 


A i A. 

18 9 18 


Table 3.1: A class of symplectic RK methods of order 4 retrieved by Gaussian quadrature. 


Theorem 3.2. If a csPRK method (12.611 with coefficients {Ar^a, Br, Cr] B^-, Cr) satisfying 

Br =Br = 1, 

Cr = Cr = T, 

A — ^ — A 

■^T, (J - 

and Ar^a is given such that both ^{r]) and (D(C) hold, then this csPRK method is symplectic and of 
order at least 

2 min(r/, (^) + 1. 

Proof. Since Br = Br ^ I,Cr = Cr = t, this implies that B{oo) holds. In such a case, note that 
Ci(r^) and (D(C) are equivalent to C{r]) and T’((^), respectively. Next, by combining Theorem 12.31 and 
Lemma l3.1l it gives the final result. □ 

One can then easily construct symplectic PRK methods by means of Theorem 13.21 and Theorem 
[23] based on numerical quadrature formulae, we provide the following corollary. 

Corollary 3.2. Assume Ar^a is a bivariate polynomial of degree d'^ in r and degree d^ in a, 
Ar^a = 1 — Afj^r, o,nd the quadrature formula {bi,Ci)l^.^ is of order p. If (^r,o-, Br, Cr) satisfies 
Br ^ 1, Cr = T (then iB(oo) holds) and both €{r]), hold, then the classical PRK method with 

coefficients {bjAci^cj,bi, Ci', bjAci^cj,bi, Ci) is symplectic and of order at least 

min(p, 2 q: + 1), 

where a = min(? 7 , C^, p — dl, p — d^). 

Proof. The symplecticity-preserving property of the associated PRK method can be easily seen 
from the following fact 

bfibjAc^^cj) + bj{biAcj^ci) — bibj, i,j = 1, • • • ,r. 

Note that Ar^ o- = 1 — r is a bivariate polynomial of degree d^ in r and degree d'^ in a, thus 
inserting dP = d^ and = dP into the formula given in Theorem 12.51 gives the final result. □ 


Example 3.2. Let Br = Br ^ 1, Cr = Cr = t and take 

s—1 s—2 


AT,a = ^ + '^^,+iP,+i{T)Pfia) - '^f,+iP,+i{a)PfiT), 


L=0 


L=0 


s-2 


s-1 


(3.8) 


Ar,a- = 1 - A^^r = - + ^^i+iPt +1 (r)P,(cj) - ^^,+iR,+i(cr)P,(r), 


t=0 


L=0 
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5-v^ 

1 5 

2 

2 VT 5 

5 2^l5 

5-715 

5 1 715 

36 90 

2 

7l5 

5 

2715 

10 


9 

45 

36 45 


9 

15 

36 

45 

1 

5 VIE 

36 24 


2 

5 VTE 

1 

5 1 715 

36 36 


2 

5 

715 

2 


9 

36 24 

2 


9 

36 

36 

5+v^ 

5 1 2 VT 5 

36 45 

2 

1 2 VT 5 
45 

5 1 715 

36 90 

5+715 

5 1 2715 

36 45 

2 

I 715 
15 

5 

715 

10 

9 

10 

9 

36 

90 


5 


4 

5 


5 


4 


5 


18 


9 

18 


18 


9 


18 

Table 3.2: A new 

symplectic PRK method of order 4 retrieved by Gaussian quadrature. 



0 

6-76 

10 

0 

19 76 

150 225 

0 

1 , 76 

4 72 

0 

67 19776 

300 1800 

0 

6-76 

10 

1 

9 

1 

9 

1 I 76 

18 72 

7 , 76 

36 72 

1 

18 

53 

180 

76 

72 

4176 

360 

6+76 

10 

19 1 76 

150 225 

67 1 19776 

300 1800 

1 76 

4 72 

6+76 

10 

1 

9 

53 1 4176 

180 360 

7 

36 

76 

72 


1 

9 

16+76 

36 

16-76 

36 


1 

9 

16+76 

36 

16 

=76 

36 


Table 3.3: A new symplectic PRK method of order 4 retrieved by Radau-Left quadrature. 


4-76 

10 

14-76 

72 

398-14776 

1800 

-7-276 

450 

4-76 

10 

1 76 

4 72 

3 3176 

20 360 

0 

4+76 

10 

398+147 76 
1800 

14+76 

72 

-7+276 

450 

4+76 

10 

3 1 3176 

20 = 360 

1 1 76 

4 72 

0 

1 

16-76 

36 

16+76 

36 

1 

9 

1 

1 76 

2 72 

1 I 76 

2 72 

0 


16-76 

36 

16+76 

36 

1 

9 


16-76 

36 

16+76 

36 

1 

9 


Table 3.4: A new symplectic PRK method of order 4 retrieved by Radau-Right qnadratnre. 


0 

0 

0 

0 

0 

0 

1 

12 

- 1+75 

24 

-1-75 

24 

5-75 

11-75 

25+75 

25-1175 

-1-75 

5-75 

1 

25-75 

25-1175 

10 

120 

120 

120 

120 

10 

12 

120 

120 

5+75 

11+75 

25+1175 

25-75 

-1+75 

5+75 

1 

25+1175 

25+75 

10 

120 

120 

120 

120 

10 

12 

120 

120 

1 

1 

5 

5 

1 

1 

1 

11+75 

11-75 

12 

12 

12 

12 

12 

24 

24 
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1 

5 

5 


1 12 

12 

12 

12 


12 

12 

12 


Table 3.5: A new symplectic PRK method of order 4 retrieved by Lobatto quadrature. 


5-715 

2 

2 2715 

1 

715 

5-715 

1 

2 

2715 

2 

715 

10 

9 

9 45 

18 

18 

10 

18 

9 

45 

9 

18 

1 

5 7l5 

2 

5 

7l5 

1 

5 7l5 


2 

5 

7l5 

2 

36 36 

9 

36 

36 

2 

36 36 


9 

36 

36 

5+715 

1 I 715 

2 1 2715 


2 

5+715 

2 1 715 

2 

_l_ 2715 


1 

10 

18 18 

9 = 45 


9 

10 

9 = 18 

9 

= 45 


18 


5 

4 


5 




4 


5 


18 

9 


18 


1 18 


9 


18 


Table 3.6: A new symplectic PRK method of order 3 retrieved by Ganssian quadrature. 
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by using some suitable quadrature formulae with s quadrature points, we can retrieve a lot of classical 
high order symplectic PRK methods including s-stage Lobatto IIIA-IIIB, Radau lA-IA, Radau IIA- 
IlA and Gauss-Legendre schemes, and the orders are 2s — 2,2s — 1,2s — 1,2s respectively. 

If we use more quadrature points, then we can get lots of new symplectic PRK schemes which 
are not found in the literature. For instance, consider s = 2, it gives the symplectic PRK schemes 
shown in Table WfMS.fA 

Example 3.3. Let Br = B^ = 1, Cr = Cr = t and take 


d j. ^ 

= - + ^^i+i(Pi+i(r)P,(cj) - P,+i(cj)P,(r)) + + 

1 1 

Ar,a = 1 -^a,r = - + ^ ^^+1 (P,+ l (r)P, (cj) - P,+ l (cj)P, (t)) - 2(2g + ^ 


(3.9) 


by using suitable quadrature formulae with s + 1 quadrature points, we can retrieve a lot of classical 
high order symplectic PRK methods including {s + l)-stage Lobatto IIIC-IIlC, Radau lA-IA, Radau 
IIA-IIA and Gauss lA-lA scheme^, and the orders are 2s, 2s + 1,2s + 1,2s + 1 respectively. 

If we use more quadrature formulae with different number of quadrature points, then we can 
get more new symplectic PRK schemes. Table \S.f^ provides a symplectic PRK scheme based on 
Gaussian quadrature for the case s = 1. 


4. Conclusions 

This paper investigates the construction of symplectic (P)RK type methods based on continuous- 
stage (P)RK methods. In the process of construction, we first discuss the general continuous-stage 
(P)RK methods with the help of the simplifying assumptions and the expansion technique of or¬ 
thogonal polynomials (i.e., Legendre polynomials), and then establish (P)RK methods in classical 
sense by means of numerical quadrature formulae. The newly obtained results are then directly 
applied to construct a special class of (P)RK methods which have symplecticity-preserving prop¬ 
erty for solving Hamiltonian systems. The line of construction of numerical algorithms combines 
the continuous and discrete mathematical theory, which forms a new and simple way for numerical 
integration of ODEs especially for Hamiltonian systems with a symplectic structure. 
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